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ABSTRACT 

We present O'.'S resolution observations of the disk around GM Aurigae with the Submilhmeter 
Array (SMA) at a wavelength of 860 /im and with the Plateau de Bure Interferometer at a wavelength 
of 1.3 mm. These observations probe the distribution of disk material on spatial scales commensurate 
with the size of the inner hole predicted by models of the spectral energy distribution. The data 
clearly indicate a sharp decrease in millimeter optical depth at the disk center, consistent with a 
deficit of ma t erial a t distances less than ~20 AU from the star. We refine the accretion disk model of 
ICalvet et al.l ()2005f ) based on the unresolved spectral energy distribution (SED) and demonstrate that 
it reproduces well the spatially resolved millimeter continuum data at both available wavelengths. 
We also present complementary SMA observations of CO J=3— 2 and J=2— 1 emission from the disk 
at 2" resolution. The observed CO morphology is consistent with the continuum model prediction, 
with two significant deviations: (1) the emission displays a larger CO J=3— 2/J=2— 1 line ratio than 
predicted, which may indicate additional heating of gas in the upper disk layers; and (2) the position 
angle of the kinematic rotation pattern differs by 11° ± 2° from that measured at smaller scales from 
the dust continuum, which may indicate the presence of a warp. We note that photoevaporation, grain 
growth, and binarity are unlikely mechanisms for inducing the observed sharp decrease in opacity or 
surface density at the disk center. The inner hole plausibly results from the dynamical influence of a 
planet on the disk material. Warping induced by a planet could also potentially explain the difference 
in position angle between the continuum and CO data sets. 

Subject headings: circumstellar matter — planetary systems: protoplanetary disks — stars: individual 
(GM Aurigae) 
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1. INTRODUCTION 

Understanding of the planet formation process is inti- 
mately tied to knowledge of the structure and evolution 
of protoplanetary disks. Of particular importance is how 
and when in the lifetime of the disk its constituent ma- 
terial is cleared, which provides clues to how and when 
planets may be assembled. While observations suggest 
that the inner an d outer dust disk disperse nearly si- 
multa neouslv fe.g. ISkrut skie et al.lll990l:IWolk fc Waited 
IT99I [Andrews fc Willia ms 200^ it is not clear which 
physical mechanism(s) drives this process, or the details 
of how it progresses. Possible dispersal mechanisms, 
of which several may come into play over the lifetime 
of a disk, in clude a drop in dust opacity due to grain 
grow th (e.g. IStrom et al.l 119891 : iDullemond fc DominikI 
l2005f ). photoevap oration of materia l by energetic stel- 
lar radiation (e.g. IClarke et al.l 20011). photophoretic ef - 
fects of gas on dust grains tera uss fc WurmI l2005f ). 
inside-out evacuation via the magnetorotational insta- 
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bility (jChiang fc Murrav-Clavl[2007[ ). and the dynamical 
interaction of giant planets with natal disk m aterial (e.g. 
iLin fc Papaloizoulll986l: iBrvden et al.lll999[ ). Observing 
the distribution of gas and dust in disks allows us to 
evaluate the roles of these disk clearing mechanisms. 

One particular cl ass of systems, those with "tran- 
sition al" disks (e.g. IStrom et al.l 119891 : ISkrutskie et al.l 
Il990( l. have become central to our understanding of disk 
clearing. These disks exhibit a spectral energy distri- 
bution (SED) morphology with a deficit in the near- 
to mid-infrared excess over the photosphere consistent 
with a depletion of warm dust near the star. The ad- 
vent of the Spitzer Space Telescope has allowed detailed 
measurement of mid-infrared spectra with unprecedented 
quality and quantity. Combined with simultaneous ad- 
vances in disk modeling th at can no w reproduce in de- 
tail the SED features (e.g. iD'Alessio et al. 1991, 120011 : 
IDullemond et all 120021 : iD'Alessio et"al..,2006 ). these ob- 
servations have revolutionized the study of disk struc- 
ture. However, such studies rely entirely on SED deficits 
whose interpretations are not unique, since effects of ge- 
omet ry and opacity can rnimic the signature of disk clear- 
ing (B oss fc Yorklfl996l : [Chiang fc Goldreichlll999D . 

Spatially resolved observations are crucial for confirm- 
ing the structures inferred from disk SEDs. High res- 
olution imaging at millimeter wavelengths is especially 
important because dust opacities are low, and the disk 
mass distribution can be determined in a straightfor- 
ward way for an assumed opacity. Millimeter observa- 
tions also avoid many of the complications present at 
shorter wavelengths, including large optical depths, spec- 
tral features, and contrast with the central star. Sev- 
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eral recent millimeter studies have resolved inner emis- 
sion cavities for disks with infrared SED de ficits through 
direct imaging observatio ns, e.g. TW Hv a (ICalvet et al.1 
2002 [ lHughes et al.l l2007f). LkHa 3 30 (Br own et al.ll2007l._ 
20081 ). and LkCa 15 f Pietu et"al.l l2007l : lE'spaillat et all 



20081 ). These observations unambiguously associate in- 



frared SED deficits with a sharp drop in millimeter opti- 
cal depth in the disk center. More information is needed 
to determine whether the low optical depth is a result of 
decreased surface density or opacity. 

GM Aurigae is a prototypical example of a star host 
to a "transitional" disk. The ^1-5 Myr old T Tauri star 
(jSimon fc Pratol[l995l iCullbring etldlfTool of spectral 
type K5 is located at a distan ce of 140 pc in the Taurus - 
Auriga molecular complex ([Bertout fc Genoval I2006D . 
and its brightness and relative isolation from intervening 
cloud material have enabled a suite of observational stud- 
ies of its disk properties. The presence of circumstellar 
du st emitting at mill i meter wavelengths was first inferred 
by iWeintraub et al.l (|1989') , and the disk structure was 
subsequently resolve d in the ^^CO J=2-l transition by 
iKoerner et al.l (|1993f ). Their arcsecond- resolution map- 
ping of the gas disk revealed gaseous material in rotation 
about the central star. Assuming a Keplerian rotation 
pattern allowed a determination of the dynamical mass 
for the central star of 0.8 Mq. Further modeling of the 
structure and dyna mics of the disk was carried out by 
iDutrev et all (|1998D . using higher-resolution ^^CO J=2~ 
1 observations. Scattered light images revealed a dust 
disk inclined by 50-56° extending to radii ^ 300 AU from 
the star (Stapclf cldt & The WFPC2 Science Team 1997; 
ISchneider et al.ll2003D . 

Efforts to model the SED of GM Aurigae have long 
indicated the presence of an inner hole, and estimates 
of its size have grown over the years as the quality of 
data and models have improved. In the early 1990s, 
the low 12 fim flux led to ^ 0.5AU estiinates of the in- 
ner d isk radius (jMarsh fc Mahorievlll992l : IKoerner et al.l 
19931). That value was later increased to 4.8 AU by 
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radiative equilibrium models, and a putative planet at a 
distance of 2.5 AU from the star was shown to be capable 
of clearing an inner hole of this extent using sim ulations 
of the relevant hydrodynamics (|Rice et al.l 20031). With 
the aid of a ground-based mid-IR spectrum. iBergin et al.l 
(|2004f ) incre ased the gap size e stimate to 6.5 AU, and 
subsequently ICalvet et al.l (p005) inferred an inner hole 
radius of 24 AU using a Spitzer IRS spectrum in com- 
binatio n with sophis t icated disk structure models. Re- 
cently, IDutrev et all (j2008f ) have argued for a 19±4AU 
inner hole in the gas distribution, using combined obser- 
vations of several different molecular line tracers. Like 
the SED-based measurements, their method is indirect: 
they use a model of the disk in Keplerian rotation to as- 
sociate a lack of high- velocity molecular gas with a deficit 
of material in the inner disk. 

We present interferometric observations at 860 fim 
from the Submillimeter Array^ and 1.3 mm from the 
Plateau dc Burc Interferometer* that probe disk material 

''' The Submillimeter Array is a joint project between the Smith- 
sonian Astrophysical Observatory and the Academica Sinica Insti- 
tute of Astronomy and Astrophysics and is funded by the Smith- 
sonian Institution and the Academica Sinica. 

Based on observations carried out with the IRAM Plateau 



on scales commensurate with the 24 AU inner disk radius 
inferred from the SED. These data allow us to directly 
resolve the inner hole in the GM Aur disk for the first 
time. We describe the observations in f|2| and present the 
dual-wavelength continuum data in We also present 
observations of the molecular gas disk in the CO J=3-2 
and J=2-l lines that allow us to study disk kinematics 
in ^3.21 We use these data to investigate d isk structure 
in the context of the SED-based models of ICalvet et al.l 
(|2005f ). described in 2) Implications for the disk struc- 
ture and evolutionary status are discussed in ^ 

2. OBSERVATIONS AND DATA REDUCTION 

The GM Aur disk was observed with the 8-element 
(each with a 6 m diameter) Submillimeter Array (SMA; 
iHo et al.l[200l in the very extended (68-509 m baselines) 
and compact (16-70 m baselines) configurations on 2005 
November 5 and 26, respectively. Observing conditions 
on both nights were excellent, with '^Imm of precip- 
itable water vapor and good phase stability. Double 
sideband receivers were tuned to a central frequency of 
349.935GHz (857 ^m), with each 2GHz-wide sideband 
centered ±5 GHz from that value. The SMA correlator 
was configured to observe the CO J=3— 2 (345.796 GHz) 
and HCN J=4-3 (354.505 GHz) transitions with a veloc- 
ity resolution of 0.18 km s~^. No HCN was detected, with 
a 3(7 upper limit of 0.9 Jy beam~^ in the 2'.'2x 1'.'9 synthe- 
sized beam. The observing sequence alternated between 
GM Aur and the two gain calibrators 3C 84 and 3C 111. 
The data were edited and calibrated using the MIR soft- 
ware package.^ The passband response was calibrated 
using observations of Saturn (compact configuration) or 
the bright quasars 3C 273 and 3C 454.3 (very extended 
configuration). The amplitude scale was determined by 
bootstrapping observations of Uranus and these bright 
quasars, and is expected to be accurate at the ~10% 
level. Antenna-based gain calibration was conducted us- 
ing 3C 111, while the 3C 84 observations were used to 
check on the quality of the phase transfer. We infer that 
the "seeing" induced on the very extended observations 
by phase noise and small baseline errors is small, <0'.'l. 
Wideband continuum channels from both sidebands and 
configurations were combined. The derived 870 /im flux 
of GM Aur is 640 ± 60 mJy. 

Additional SMA observations in the extended (28- 
226 m) and sub-compact (6-69 m baselines) configura- 
tions were conducted on 2006 December 10 and 2007 
September 14, respectively, with a central frequency of 
224.702GHz (1335 ^m). While the sub-compact ob- 
servations were conducted in typical weather condi- 
tions for this band (2.5mm of water vapor), the ex- 
tended data were obtained in better conditions similar 
to those for the higher frequency observations described 
above. The correlator was configured to simultaneously 
cover the J=2-l transitions of CO (230.538 GHz), i^CO 
(220.399GHz), and C^^^O (219.560GHz) with a velocity 
resolution of ~0.28km s~^. The calibrations were per- 
formed as above. 

GM Aurigae was also observed with the 6-element 
(each with a 15 m diameter) Plateau de Bure Interfer- 

de Bure Interferometer. IRAM is supported by INSU/CNRS 

(France ), MFC (Germany) and ION (Spain). 

^ See jhttp : //cf a-w«w. harvard. edu/$\slm$cql/mlrcook. html [ 
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Fig. 1. — Continuum emission from the disk around GM Aur at wavelengths of 860 /^m observed with the SMA (top) and 1.3mm 
observed with PdBI {bottom). The data are displayed in both the image (a) and Fourier (b) domains. In the image domain (a), the 
observed brightness distribution at each wavelength {left) is compared with the model prediction {center; see t|4.1l for model details), and 
the residuals are also shown {right). In the data and model frames, the contours are [3,6,9, ...]x the rms noise (3.5 mjy beam~^ at 860 /^m 
and 0.75 mjy beam" at 1.3mm). In the residual frame, the contours start at 2a and are never greater than 3cr. The synthesized beam 
sizes and orientations for the two maps are, respectively, 0''30 X 0''24 at a position angle of 34° and 0''43 X 0''30 at a position angle of 
35°. Two sets of axes are shown: the dotted line indicates t he po sition angle of the double-peaked continuum emission, while the solid line 
indicates the best-fit position angle of the CO emission (see i|3.2l for details). In the Fourier domain (6), the visibilities are averaged in bins 
of deprojected u-v distance from the disk center, and compared with the model prediction (red line). The inner hole in the GM Aur disk 
is clearly observed at both wavelengths, as a double-peaked emission structure in the image domain or as a null in the visibility function 
in the Fourier domain. 



ometer (PdBI) in the A configuration (up to 750 m base- 
lines) on 2006 January 15. Observing conditions were ex- 
cellent, with atmospheric phase noise generating a seeing 
disk of <0.2". The PdBI dual-receiver system was set 
to observe the 110.201GHz (2.7 mm) and 230.538 GHz 
(1.3 mm) continuum simultaneously. As with the SMA 
data, observations alternated between GM Aur and two 
gain calibrators, 3C 111 and J0528-I-134. The data were 
edited and calibrated using the GILDAS package (jPetvl 
|2005[) . The passband responses and amplitude scales 
were calibrated with observations of 3C 454.3 and MWC 
349, respectively. The derived 1.3 and 2.7 mm fluxes of 
GM Aur are 180 ± 20 and 21 ± 2 mJy. 

The standard tasks of Fourier inverting the visibilities, 
deconvolution with the CLEAN algorithm, and restora- 
tion with a synthesized beam were conducted with the 
MIRIAD software package. A high spatial resolution im- 
age of the 860 /zm continuum emission from the SMA 
data was created with a Briggs robust — 1.0 weight- 
ing scheme for the visibilities, excluding projected base- 
lines < 70 kA, resulting in a synthesized beam FWHM of 
O'.'SO X 0'.'24 at a position angle of 34°. A similar image of 
the 1.3 mm continuum emission with a synthesized beam 
FWHM of 0'.'43 X 0'.'30 at a position angle of 35° was 
generated from the PdBI data using natural weighting 
(robust — 2.0). Tabic [T] summarizes the line and contin- 
uum observational parameters. 

3. RESULTS 

3.1. Millimeter Continuum Emission 

Figure [1] shows the results of the SMA and PdBI con- 
tinuum observations in both the image and Fourier do- 
mains. The presence of an inner hole in the GM Aur disk, 
as predicted by models of the SED, is clearly indicated 
both by the double-peaked emission structure in the im- 
age and by the null in the visibility data. The double- 
peaked emission structure points to a deficit of flux near 



the disk center; the null in the visibility function, or the 
location at which the real part of the visibilities change 
sign, similarly reflects a decrease in flux at small angular 
scales. The resolution of the 2.7 mm data from the PdBI 
was insufficient to provide about the inner hole. 

The maps in the left panel of Fig. [1] show a double- 
peaked brightness distribution at both wavelengths, with 
peak flux densities of 59 ± 4 mJy beam~^ at 860 /Ltm and 
16.6 ± 0.3 mJy beam~^ at 1.3 mm. For all but t he most 
edge-on viewing geometries (e.g. lWolf et a l. 2008). a con- 
tinuous density distribution ext ending in to the d ust de- 
struction radius (-0.05-0.1 AU; llsella et al.ll2006( ) would 
be expected to result in a centrally-peaked brightness dis- 
tribution. In the case of GM Aurigae, the double-peaked 
emission structure is a geometric effect due to the trun- 
cation of disk material at a much larg er radius, viewed at 
an in termediate inclination of 50-56° (iDutrev et al.llT998l . 
|2008[ ) : the region of highest density is near the inner disk 
edge, with a large column density of optically thin ma- 
terial in this ring effectively generating limb brightening 
at the inner edge of the outer disk, at two points along 
the disk major axis. 

The size of the inner hole can be roughly estimated 
by the separation of the emission peaks, although the 
peak separation will also depend on the brightness of the 
directly-illuminated inner edge of t he outer disk r elative 
to the extended disk component (jHughes et al.l 12007) . 
The separation of the peaks in the 860 /xm image is 0'.'38± 
0'.'03, corresponding to a physical diameter of 53 ± 4 AU 
(radius 27 ± 2 AU) at a distance of 140pc. A position 
angle of 66° is estimated by the orientation of a line that 
bisects the two peaks, although a more robust value of 
64° ± 2° is derived in §4.11 below. Since the peaks are 
not distinctly separated in the 1.3 mm image, the same 
estimate cannot be made, but the position angle is clearly 
consistent with that derived from the 860 fim visibilities 
and indicated by the perpendicular dashed lines in Fig. m 
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TABLE 1 

Observational parameters for GM Aur 



Parameter 


"CO J=3-2 


i2cO J=2-l 


860 Atm 


Continuum 
1.3 mm 


2.7 mm 


Rest Frequency (GHz) 


345.796 


230.538 


349.935 


230.538 


110.201 


Channel Width 


0.18kms-i 


0.28kms-i 


2x2 GHz 


2 X 548 MHz 


548 MHz 


Beam Size (FWHM) 


2'.' 2x1'.' 9 


2"lxl^'4 


0y30x0'.'24 


0'.' 43x0'.' 30 


0^' 93x0'.' 60 


PA 


14° 


56° 


34° 


35° 


31° 


RMS noise (mjy beam"-"^) 


310 


90 


3.5 


0.75 


0.25 


Peak Flux Density (mjybcam"^) 


6700 ± 300 


2400 ± 100 


59 ±4 


16.6 ± 0.8 


10.3 ±0.3 


Integrated Continuum Flux (mjy) 






640 ± 60 


180 ± 20 


21 ± 2 


Integrated Line Intensity (Jykms"'^) 


29 


37 









The presence of an inner hole is also evident from the 
visibilities displayed in the right panel of Fig. [1] The 
real part of the complex visibilities have been averaged 
in concentric annuli of deprojected (u, v) distance from 
the disk center. For de tails of the deprojection pro- 
cess, see L av et all (119971) . As discussed in the appendix 
of iHughes et al.1 (|2007f l. the presence of a null in the 
visibility function indicates a sharp decrease in flux at 
a radius corresponding roughly to the angular scale of 
the null position. The precise position of the null de- 
pends primarily on the angular size of the inner hole, 
but also on the radial gradients of the surface density 
and temperature distribution and the relative brightness 
of the directly illuminated wall at the inner edge of the 
outer disk. In a standard power-law parameterization, 
the disk temperature T and surface density S vary in- 
versely with radius as S (x and T c>c Ne- 
glecting the emission from the wall and assuming stan- 
dard values oi p = 1.0 and q = 0.5, expected for a 
typic al viscous disk with constant a (jHartmann et all 
Il998f l and consistent with p revious studies of the GM 
Aur disk jP utrey et al.lll998l : [Andrews fc WiUiams 2007; 
IHughes et al. 2008Y we may obtain a rough estimate of 
the size of the inner hole using the obser ved null posi- 
tion and Eq. A9 from IHughes et all (|2007D : 7^„ull(kA) = 
(1 AU/i?hoic)(£'sourcc/100 pc)[2618 + 1059(p + q)]. A 
polynomial curve fit to the visibilities yields a null posi- 
tion of 190 kA at 860 /j,m and 224 kA at 1.3 mm, which 
correspond to inner hole radii of 31 and 26 AU, respec- 
tively. However, these estimates are uncertain to within 
~30%, as the data are consistent with a broad range of 
null positions. We therefore turn to a more sophisticated 
modeling procedure described in i )4.1l below. 



3.2. CO Channel and Moment Maps 

Figures I2H11 display the new SMA observations of CO 
emission from the GM Aur disk. Figures [2] and |3| show 
channel maps with contours starting at twice the rms 
noise level and increasing by factors of V2, while Figure |4| 
displays the zeroeth (contours) and first (color) moments 
of the data: these are the velocity-integrated intensity 
and intensity-weighted velocities, respectively. The peak 
flux density is 6.7±0.3 Jybeam"^ in the CO J=3-2 line 
and 2.4±0.1Jybeam~^ in the CO J=2-l line, with inte- 
grated fiuxes of 9.4 Jykms"^ and 21.2 Jykms"^, respec- 
tively (although emission from extended ambient cloud 
material is likely to increase the CO J=2-l integrated 
flux over that originating from the disk alone) . The chan- 
nel and moment maps are broadly consistent with the 
expected kinematic pattern for material in Keplerian ro- 
tation about the central star, substantially inclined to 



our line of sight (as in iDutrev et al.l Il998t I Simon et al.l 
I2OOQ ). 

The short-baseline spatial frequencies in the (u, v) 
plane provided by the subcompact configuration of the 
SMA during our observations of the J=2— 1 transition 
are sensitive to emission on the largest spatial scales. 
These short antenna spacings reveal the severity of the 
cloud contamination to an extent not possible with pre- 
vious data. The contamination is evident as an ex- 
tended halo around the disk emission in the central chan- 
nels of the ,7=2—1 channel maps near LSR velocities 
of 5-6 km s^^ (Fig- El)- It is also evident in the mo- 
ment map (Fig. [4]) as an elongation of emission near 
the systemic velocity (green-yellow) to the northwest 
along the disk minor axis. This contamination indi- 
cates that caution must be exercised when deriving kine- 
matic information from the CO lines, particularly the 
central channels. Spatial filtering by the interferometer 
does not ameliorate cloud contamination in an abundant, 
easily-excited, high-optical depth tracer like C0J=2-1. 
The J=3— 2 line appears less contaminated than J=2— 1 
(Figs, m and ID), although similarly short antenna spac- 
ings (8-43 m) are not present in this data set. Neverthe- 
less, we expect less cloud contamination in the J=3— 2 
transition, since the temperature of the cloud will be 
lower than that of the disk and will therefore populate 
the upper rotational levels of the CO molecule less effi- 
ciently. The cloud contamination prevents detection of 
self-absorption in the central channels of the CO J=2- 
1 channel maps along the near (northwest) edge of the 
disk (as deterniined b y scattered light obse rvations; see 
ISchneider eraI]|2003D . iDutrev et all (|1998D report self- 
absorption along the southeast edge, but our observa- 
tions suggest that this brightness asymmetry may be due 
to cloud contamination. It is also possible that the con- 
tamination is due to a residual envelope, although we 
are unable to determine the large-scale structure of the 
extended line emission with our interferometric data. 

In all figures, the disk orientation based on the posi- 
tion angle of 64° derived from the continuum emission 
(Fig. [1] and M.l^ is plotted over the CO emission as a 
set of crossed dashed lines, with the relative extent of 
the major and minor axes (based on the inclination an- 
gle of 55°) indicated by the length of the perpendicular 
lines. The position angle of 51° derived by IDutrev et al.l 
()1998l ) from fitting the CO J=2— 1 emission, consistent 
with our own J=3— 2 and J=2— 1 observations, is illus- 
trated by the solid line. Note that the position angle of 
the CO emission difl:ers slightly from the position angle 
of the continuum emission, by 11° ± 2° (see ^TA^. The 
trend is clear for both transitions, but more obvious in 
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Fig. 2. — Channel maps of CO J=3— 2 emission from the GM Aur disk. Contour levels start at 0.61 Jy (2 times the rms noise) and 
increase by factors of \/2. LSR velocity is indicated by color and quoted in the upper right of each panel. The synthesized beam (2'.'2xl'.'9 
at a PA of 14°) and physical scale are indicated in the lower left panel. Two sets of axes are shown: the dotted line indicates the position 
angle of the double-peaked continuum emission, while the solid line indicates the best-fit position angle of the CO emission. 




Fig. 3. — Channel maps of CO J=2— 1 emission from the GM Aur disk. Contour levels start at 0.17 Jy (2 times the rms noise) and 
increase by factors of \/2. LSR velocity is indicated by color and quoted in the upper right of each panel. The synthesized beam (2" 1 X 1''4 at 
a PA of 56°) and physical scale are indicated in the lower left panel. Two sets of axes are shown: the dotted line indicates the position angle 
of the double-peaked continuum emission, while the solid line indicates the best-fit position angle of the CO emission. Cloud contamination 
is evident in at least the central four channels. 



the less-contaminated J— 3— 2 transition. Note that the 
position angle for the CO emission is derived entirely 
from the rotation pattern (evident in the isovelocity con- 
tours) and not fr'om the geometry of the integrated CO 
emission: the integrated emission appears to match the 
position angle from the continuum emission reasonably 
well. We do not observe th e isophote twi sting in inte- 
grated CO emission seen bv iDutrev et al.l ([1998). The 
cloud contamination and differences in antenna spacings 
may play a role. 

4. DISK STRUCTURE MODELS 

4.1. Updated SED Model 

Here we revisit the broadband SED m odeling of 
GM Aur presented by ICalvet et all (|2005f ). Taking 
into consideration new observational constraints at sub- 
millimeter and millimeter wavelengths, we use the irra- 
diate d accretion disk models of D'Ales sTo et al.l (|2005l 
[200l) to re-derive the properties of the outer disk of GM 
Aur and its inner, truncated edge or "wall." Our grain- 
size distribution follows a power-law of a~^'^ , where 



a is the grain radius. We assume ISM-sized grains 
in the upper layers of the disk an d accordingly adop t 
amiri=0.005 /im and a„iax—^-2^ /^m ( Draine fc Leell984) . 
Closer to the disk midplane grains have a maximum 
size of 1 mm. Input parameters for the outer disk in- 
clude the stellar properties, the mass accretion rate, the 
viscosity parameter (a), and the settling parameter (e) 
which measures the dust-to-gas mass ratio in the up- 
per layers of the disk relativ e to the standard dust- 
to-gas mass ratio. Following ICalvet et al.l (|2005f ). we 
adopt the same extinction, distance, inclination, dust 
grain opacities, and stellar properties (i.e. luminosity, 
radius, and temperature; see Table [2]). We use a mass 
accretion rate of 7.2x 10~^ M p yr^-'^ which was derive d 
using HST STIS spectra bv llnglebv fc CalvetJ jlOOi), 
in contrast to the value of 10~^Mf7^y r ~^ de rived from 
veiling measurements in ICalvet et al.l (|2005D . We as- 
sume an outer disk radius of 300 AU, which matches 
the observed extent of scattered light from the dust disk 
(jSchneider et al.l I2003D and previo us fits to the contin- 
uum emission ( Hughes et al.ll2008h . as well as the short- 



6 



Hughes et al. 



<i 



-1 — 1 — 1 — I— T- 



GM Aur 
CO 3-2 





200 AU 



2 
A« (") 



GM Aur 
CO 2-1 




200 AU 



Fig. 4. — Zeroeth (contours) and first (colors) moment map of 
the CO J=3-2 (top) and J=2-l (bottom) data in Figs. [2] and |3] 
The dotted line indicates the position angle of the double-peaked 
continuum emission, while the solid line indicates the best-fit po- 
sition angle of the CO emission. The zeroeth moment contours 
are well aligned with the latter, while the isovelocity contours of 
the first moment map are more consistent with the former. Cloud 
contamination is evident in the CO J=2— 1 map in the northwest 
region along the disk minor axis. 

baseline data presented here. 

In order to reproduce the outer disk component of 
the SEP, we vary e and a (Figure [5]). As described in 
(jCalvet et al.l I2005D , a effectively determines the mass 
surface density distribution and therefore the disk mass, 
which is best reflected by the longest-wavelength SED 
points. The value of e has the greatest effect on the 
slope of the SED beyond 100 /im. With the new mil- 
limeter data we find e=0.5, indicating less settling than 
reported previously. We also find a=0.002 and a more 
massive outer disk of O.I6M0. This mass is significantly 
larger than an estimate based on the 860 /xm and 1.3 mm 
flux m easurements using opacities from iBeckwith et al.l 
(|1990l ). which yields -O.O4M0, and is only marginally 
Toomre stable at 300 AU (Q'-l.l). The outer disk model 
uses a n opacity of ^0.1 cm^ g^^ at 1 mm (jP'Alessio et alJ 
1200 It) which i s about four ti mes lower than that de- 
rived from the IBeckwith et al] (|T990) opacities, account- 
ing for the discrepancy in mass. Within the inner disk 
hole, there are 1.1x10^^"'^ Mq of optically thin small dust 
grains, which account for the 10 /im emission and the 



near-IR excess. The mass in solids could be much larger 
than this mass if pebbles, rocks, or even planetesimals 
have grown in the inner disk, since they woul d have a neg- 
ligible opacity in the near-IR. We note that ICalvet et al.l 
(|2005') reports the mass of the dust as 7x 10"^° Mq; this 
is actually the mass of the gas within the hole, assuming 
the standard dust to gas mass ratio. The gas mass could 
be signiflcantly larger, depending on the total amount 
of solids and the actual ratio, but these are poorly con- 
strained by existing data. 

We vary the temperature of the wall to best repro- 
duce the data. The radius of the wall is set by the tem- 
perature and dust composition, and the wall's height 
is set by the disk scale height. We assume that the 
wall is axisymmetric and composed of relatively small 
grains, as well as vertically flat in order to reproduce 
the rapid rise of the mid-IR excess at wavelengths be- 
yond 10 iim. We ado pt th e dust compo s ition used in 
iD'Alessio et"an (|2005[ ) and ICalvet et al.] ()2005f ). The 
maximum grain size is adjusted from ISM sizes to re- 
produce the shape of the IRS spectrum as necessary. At 
short wavelengths, larger grains have smaller opacities 
than ISM-sized grains. Therefore, at a given tempera- 
ture large grains will be at smaller radii than ISM -sized 
grains as per Eqn. 12 of I D'Alessio et al.l (|2005[ ). The 
derived size of the inner hole varies somewhat depend- 
ing on whether the SED or the resolved millimeter vis- 
ibilities are included. Fitting only the broadband SED 
and neglecting the resolved millimeter-wavelength data, 
the wall is located at 26 AU and has a temperature of 
130K and a height of ^2 AU with maximum grain size 
amax=0-'25 ^m. (Fig.[5l left p anel) . The ra dius of the wall 
differs by -2 AU from .Calvet et all (|2005l ). since here we 
take Lace ~ GMM/R assuming magnetospheric accre- 
tion while ICalvet et all (2005) uses Lace ^ GMM/2R as 
per the boundary layer model. We also adopt a different 
mass accretion rate. 

In order to compare the SED model with the resolved 
continuum data, it is necessary to fix the disk geometry. 
As listed in Table [21 we adopt an inclination of 55°, i n 
order to maintain consistency with ICalvet et ahl (|2005| ). 
However, the position angle is poorly reproduced by the 
value of 53.4° ± 0.9 ° that is the w eighted average of fits 
to the CO emission (iDutrev et al . 1998, 2008, see Fig.[l]). 
To derive a more appropriate position angle, we generate 
a sky-projected image from the disk model and use the 
MIRIAD task uvmodel to sample the image at the same 
spatial frequencies as the data. We compare these model 
visibilities with the observed 860 /im visibilities (which 
have the finest resolution). We repeat this process for a 
range of position angles and calculate a value com- 
paring each set of model visibilities with the data. Using 
this method, we fit a position angle of 64° ± 2°, which 
differs by 1 1° zfc 2° from the pos i tion a ngle of the CO disk 
derived bv Du trev eFall ([1991, [200l) . 

When considering the resolved millimeter-wavelength 
visibilities, a disk with a 20 AU hole reproduces the emis- 
sion much better (Fig. El right panel, and Fig. [H center 
panels). Using the same comparison of visibilities as 
described in the previous paragraph, the 20 AU model 
represents a Str improvement over the 26 AU model, 
which significantly underpredicts the amount of flux pro- 
duced close to the star. This 20 AU hole has a wall with 
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a temperature of 120 K, a height of 1.4 AU, and maxi- 
mum grain size amax=5 fim. For neither the 20 AU nor 
the 26 AU model does the wall contribute significant con- 
tinuum emission at the wavelengths and spatial scales 
probed by our data. The main discrepancy between the 
fits to the SED and the millimeter visibilities occurs be- 
tween wavelengths of ~20-40/im where the 20 AU hole 
model overpredicts the flux. However, the SED morphol- 
ogy in this region is likely sensitive to the properties of 
the wall at the inner disk edge, which are not well known 
and are not constrained by our data. It is also possible 
that the composition of the grains, particularly whether 
the silicate and graphite form composite grains or are 
separated, can affect the temperature and therefore the 
mid-IR morphology of the wa ll component of the SED 
p'Alessio et al.ll2009. in prepD . Since our focus is on the 
interferometric millimeter- wavelength data, we adopt the 
model with a 20 AU inner hole for the remainder of the 
analysis. Figure [1] compares this model with the data 
in the image plane (center panel) and in the visibility 
domain (red line in the right panel). The agreement is 
excellent, and the residuals are less than 3ct within the 
2" box shown. 

The flux density of the eastern peak of the 860 im- 
age is 50mJybeam^^, while that of the western peak is 
59mJybeam^^. The corresponding peaks in the model 
images are 49 and 50mJybeam~^, respectively. Given 
the rms noise of 3.5 mJy beam~^ , these values are consis- 
tent with no flux difference and hence axially symmetric 
emission from the inner disk edge. The positional ac- 
curacy of the data and knowledge of the stellar proper 
motion are insufflcient to determine whether or not the 
emission peaks are equally offset from the star. This re- 
sult may be contrasted with the strong asymmetries ob- 
served bv lBrown et al.l (|200^ in their observations of the 
inner hole in LkHa 330, although these data are missing 
short antenna spacings present in the GM Aur data that 
may dilute asymmetries. However, as in the case of LkHa 
330, we find that the GM Aur continuum presents a sharp 
contrast in brightness between the inner and outer disk, 
reflected by the null in the visibility function and the 
strong agreement between the data and the model con- 
taining an inner hole. The 1.1 x 10~^^ Mq of dust within 
the central hole in the model implies a reduction in the 
mass surface density of small grains of at least 6 orders of 
magnitude at 1 AU relative to a continuous model of the 
dust disk, indicating that the data are consistent with an 
inner disk region that is essentially completely evacuated 
of small grains. 

4.2. Comparison with CO Observations 

In order to compare the gas and dust properties of the 
GM Aur disk, we used the SED-based model described 
above to generate predicted CO J=3— 2 and J=2— 1 
emission. We assume that gas and dust are well mixed, 
with a uniform gas-to-dust mass ratio of 100 (neglect- 
ing the complication of dust settling) and a constant CO 
abundance relative to H2 of 10^^, which is required to 
reproduce the peak CO J=2— 1 flux. We also add mi- 
croturbulence with a FWHM of 0.17kms~ ^ throughout 
the outer disk, as derived bv lDutrev et al.l ([1998). This 
is comparable to the 0.18 kms~^ spectral resolution of 
the data and does not affect our determination of the 
disk geometry. Due to the position angle differences ev- 



TABLE 2 
Stellar and Model Properties 





1.1 


Ft ( Rr^) 


1.5 


l^j 


^ / oU 


M (M0yr-i) 


7.9 X 10-3 


Distance (pc) 


140 


Av 


1.2 


Inclination (dcg) 


55 


Optically Thick Wall 




R^all (AU) 


20 (26) 


amin (/im)i 


0.005 


a.max {A»ni)3 


5 (0.25) 


T^all {Kf 


120 (130) 


Z^all {AC/)3.4 


1.4 (2) 


Optically Thick Outer Disk 


Rd,out {Auy 


300 


e3 


0.5 


0? 


0.002 


Ma(MQ) 


0.16 



ident between the continuum emission in Fig. [T] and the 
central channels in F ig. [2l we also adjust the position 
angle to 51° (as in iDutrev et al.l fT998) . Finally, we 
note that with an outer radius of 300 AU, the contin- 
uum model severely underpredicts the CO emission at 
large radii, as expec ted for a model wit h a sharp cut- 
off at its outer edge (jHughes et al.ll20()8l ). We therefore 
extrapolate the model to 525 AU to match th e spatial ex- 
tent of the CO emission (jDutrev et al.ll 1998 V While this 
larger CO model no longer matches perfectly the con- 
tinuum emission for the shortest baselines, based on the 
prediction assuming a constant gas-to-dust mass ratio, it 
retains the kinematic and thermal structure of the small- 
scale continuum model. In order to consistently solve for 
the level populations and generate sky-projected images 
in the CO lines, we use the Monte Carlo radiative tr ansfer 
code RATRAN (jHogerheiide fc van der fagi2000t l. We 
then use the MIRIAD task uvmodel to sample the model 
image at identical spatial frequencies to those present in 
our interferometric CO data set. 

Figure [5] compares the predicted CO emission from the 
extended SED model (right) with the observed emission 
from the GM Aur disk (left) for the J=2-l (top) and 
J=3— 2 transitions. It is clear that the velocity pattern in 
the disk is consistent with Ke plerian rotation ( as previ- 
ously noted by IKoerner et al.lll993tlDutrev et aLl ll998). 
and that the SED-based model is capable of reproducing 
the basic morphology of the CO emission. 

The primary difference between data and model is the 
CO J=3-2/J=2-l line ratio: the disk structure model 
that reproduces the peak flux density of the J=2— 1 tran- 
sition underpredicts the peak J=3— 2 flux by 30%. This 
difference may be attributed to a ~10K difference in 
temperature between the gas and dust in the upper lay- 
ers of the GM Aur disk that are probed by these op- 
tically thick CO lines. While the vertical temperature 
gradient of the dust in the model is flxed by the SED, 
a relative increase in gas temperature would populate 
the upper rotational transition of the molecule more ef- 
flciently and produce more J=3— 2 emission relative to 
J=2— 1. The temperature and the CO abundance are 
also somewhat interdependent, since the CO abundance 
sets the vertical location, and therefore the temperature. 
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Fig. 5.— Model of the SED of GM Aur using the method of lD'Alessio et aLI pOOSl . 120061 ') . The final model of the SED alone has an inner 
disk hole of 26 AU {left), while the model that best reproduces the resolved millimeter- wavelength visibilities has a hole of radius 20 AU 
(right). See jplfor mod el details. We show optical ( open circles; [Kcnvon & Hartmann 1995), 2MASS (closed cir cles), IRAC (open squares; 
IHartmann et all 12005"). and IRAS (closed squares; Weaver & Jones 1992) data and a Spitzer IRS spectrum (Calvet et al.l 120051') . Open 
pentagons repre sent millimeter observations obtained from Andrews & Williams (2005); Beckwith & Sargent (1991); Du trev et al.l (l998); 
IKitamura et al.l 1^002); Koe rner et al.i H1993I '): iLoonev et al.i(j200 0l: Rodm ann et al.j H2006i '): iWeintraub et al.i (|1989r '). Closed pentagons are 
from this work. The final model (solid line) includes the following components: stellar photosphere (dotted line), optically thin dust region 
(long-dashed line), disk wall (short-long dashed line), outer disk (dot-dashed line). The peak at ^l/im from the wall emission is due to 
scattered light. While the 20 AU model does not fit the IRS spectrum as well between ~20-40 /xm as the 26 AU model, it reproduces the 
millimeter continuum emission very well at both 860 lira and 1.3mm (Fig.[TJ. 

of the r=l surface from which most of the hue emission 
originates. An increase in temperature would therefore 
also vary the anomalously low CO/H2 ratio necessary to 
reproduce the J— 2—1 flux. Such line ratio differences 
have been previously obs erved in the disk around TW 
Hya (iQi et al.ll20dll2006h . and may be due to additional 
heating of gas in the upper disk by such processes as x- 
ray and UV irradiation, dissocia tive or mechanical heat- 
ing (e.g.[Gia ssgold et al. 2004; K amp fc Dullemondll2004l : 
iNomura et'al.ll2007l) 

Nevertheless, while the flux levels vary between the 
data and model prediction, the similarity in morphology 
makes it clear that the overall disk structure is consis- 
tent between the molecular gas traced by CO and the 
model based on dust traced by continuum emission and 
the SED. The only other significant difference between 
the two is in the position angle of the emission, which 
differs by ~11°. The implications of this result are dis- 
cussed in ij5.2l below. 



5. DISCUSSION 

5.1. Inner Disk Clearing 

The resolved millimeter continuum observations of the 
GM Aur system are consistent with the prediction from 
the SED model. Models of the observed 860 ^m and 
1.3 mm maps in conjunction with the SED and Spitzer 
IRS spectrum, give a value of ^20 AU for the extent 
of this inner cleared region. The inference of an inner 
hole of this size from the SED and resolved millimeter 
visibilities is consistent with recent millimeter-wave ob- 
servations of rotational transitions of CO isotopologues 
from the GM Aur disk that provide spectroscopic ev- 
iden ce for a diminished density of cold CO within 20 
AU ()Dutrev et al.l [20081 ) . However, other observations 
indicate that this region cannot be entirely devoid of 



gas. ISalvk et "ahl ()2007[ ) detect CO rovibrational emis- 
sion originating from hot gas at radii near ~ 0.5 AU, 
from which they infer a total gas mass in the inner disk 
of ~ 0.3 Mg. Measurements of the Ha l inewidth imply 
an a c cretion rate of ^ 10~^ M© yr~^ ([White fc Ghej 
[20011: llnglebv fc Calvetll2009f) : accretion at this rate re- 
quires a steady supply of gas from the inner disk. The 
SED model also requires 3 x 10"'' lunar masses of dust in 
the inner disk, to account for the 10 nm silicate feat ure 
and slight near- to mid-IR excess ([Calvet et al.l[2005[ ). 

A wide variety of mechanisms has been invoked to ex- 
plain the low optical depth of the central regions of tran- 
sition disks (see e.g. [Naiita et al.l [20071 and references 
therein), each with different implications for planet for- 
mation and the process of evolution between the primor- 
dial and debris disk stages. The available measurements 
of properties of the inner hole in the GM Aur disk allow 
us to evaluate the plausibility of each mechanism as the 
driver of disk clearing in this system. 

Grain Growth - The agglomeration of dust into larger 
particles should proceed faster in central regions where 
relative velocities of particles are faster and surface den- 
sities are higher. This would produce a drop in opacities 
associated only with the inefficiency of emi ssion of large 
grains at the o bserved wavelengths (e.g. [Strom et al.l 
1989; Dullemon d fc Dominik 2 005,) . However, this pro- 
cess is inconsistent with the clearing of CO from the 
central region observed bv [Dutrev etall ([2008( 1 . as gram 
growth should proceed without diminishing the gas den- 
sity. Grain growth is also somewhat inconsistent wi th the 
steep submillimeter slope observed by Rodmann et al.l 
([20061 ) for the GM Aur system. The value inferred for 
the millimeter wavelength slope a of 3.2 is the steepest in 
their sample of ten T Tauri stars, and is typical of a grain 
population that has undergone little growth, with grain 
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Fig. 6. — Position-velocity diagram comparing the molecular line observations (left) with the predicted (right) CO J=2— 1 (top) and CO 
J=3-2 (bottom) emission from the GM Aur disk, assuming a standard gas-to-dust mass ratio of 100. The plots show the brightness as a 
function of distance along the disk major axis, assuming a position angle of 51°. Contours are [2,4,6,...] times the rms flux density in each 
map (0.17 and 0.61 Jy beam~^, respectively). The dotted line shows the expected Keplerian rotation curve for a star o f ma ss 0.84 Mq. 
The outer radius of the model has been extended to 525 AU to reproduce the extent of the molecular gas emission (see i|4.2l for details). 
The CO morphology is consistent with the SED-based model, with the exception of the line ratio: the model that best reproduces the peak 
flux of the CO J=2-l line underpredicts the CO J=3— 2 brightness by 30%. 

size Gmax ^ 1 Him. Furthermore, the original SED model 
and the submillimeter visibilities both independently in- 
dicate a sharp decrease in surface density or opacity near 
24 AU, while grain growth and dust settling are pre- 
dicted to be a continuous process and so should display 
a mo re gradual transition between the inner and outer 
disk (IWeidenschilling et al.lllQQTtlDullemond fc DominiS 
[20051) . 

Photoevaporation - Another proposed process to 
generate inside-out clearing of protoplanetary disks 
is photoevaporatio n via the "UV switch" mechanism 
(jClarke et al.ll200i[) . In this scenario, high-energy pho- 
tons from the star heat the upper disk layers, allow- 
ing material to escape the system at a rate that grad- 
ually diminishes the disk mass, while most of the disk 
mass drains onto the star via viscous accretion (e.g. 
iHartmann et al.lll998f l. Once the photoevaporation rate 
matches the accretion rate near 1 AU and prevents resup- 
ply of material from the outer disk, the inner disk will de- 
couple and drain onto the star within a viscous timescale, 
leaving an evacuated central region surrounded by a low- 
mass outer disk that will then rapidly disperse. As noted 



bv [Alexander fc Armitaed (|2007f ). the properties of the 
GM Aur system are inconsistent with a photoevaporative 
scenario because the large mass of the outer disk should 
still be sufficient to provide a substantial accretion rate 
to counteract the photoevaporative wind. Furthermore, 
the measured accretion rate is high enough that within 
the framework of the photoevaporation scenario, it would 
only be observed during the brief period of time when the 
inner disk was draining onto the star. Photoevaporation 
may yet play a role in clearing the outer disk of its re- 
maining gas and dust, but it cannot explain the current 
lack of inner disk material. 

Inside- Out MRI Clearing - The magnetorotational in- 
stability operating on the inner disk edge may also drive 
accretion and central clearing, although it should be 
noted that this is purely an evacuation mechanism: it 
can only take hold after the generation of a gap by some 
other means. Nevertheless, given the creation of a gap, 
MRI clearing is predicted to operate in systems like GM 
Aur whose outer disks a re still too massive for p hotoe- 
vaporation to dominate llChiang fc Murrav-Clavl [2007.). 
The observed depletion of CO interior to 20 AU radius 
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(jPutrev et al.|[2008l ) is consistent with this theory, which 
predicts a total gas mass depletion of order lOOOx inte- 
rior to the rim radius relative to the extrapolated value 
from the outer disk power law fit, normalizing to the to- 
tal disk mass of O.I6M0. This theory is consistent with 
the substantial accretion rate of the GM Aur system, 
yielding a value of a of 0.005, only shghtly greater than 
the de rived value of 0.002 from the model. iSalvk et al.l 
( |2007l ) estimate a gas-to-dust ratio of ~ 1000 in the inner 
disk, roughly 10 times greater than that of the outer disk, 
which is consistent with the prediction of the inside-out 
MRI evaporation scenario that flux from the star should 
promote blowout of small dust grains by radiation pres- 
sure, substantially clearing the inner disk of dust even as 
the gas continues to accrete onto the star. However, it 
is difficult to reconcile this with the substantial popula- 
tion of /Ltm-size grains that must be present in the inner 
disk to account for the 10 /xm silicate feature in the IRS 
spectrum. It is also important to consider the source of 
the requisite initial gap in the disk. 

Binarity - The dynamical influence of an unseen stel- 
lar or substcUar companion would also cause clearing of 
th e inner disk. A notable example is the recent result 
by [Ireland fc Kraiii ()2008D demonstrating that the inner 
hole in the transition disk around CoKu Tau/4 is caused 
by a previously unobserved companion. There are rela- 
tively few constraints on the multiplicity of GM Aur at 
the < 20 AU separations relevant for the inner hole. Ra- 
dial velocit y studies with km s~^ precision do not not e 
variability (jBouvier et al.l 119861 : iHartmanii et al.lll986D. 
ruling out a close massive companion. As iDutrev et al.l 
((2OO8) discuss, the stellar temperature and dynamical 
mass from the disk rotation combined with the H-band 
flux place an upper limit of ^^0.3 on the mass of a 
companion. Interferometric aperture-masking observa- 
tions with NIRC2 that take advantage of adaptive optics 
on the Keck II telescope place an upper limit of ^^40 
times the mass of Jupiter on companions with separa- 
tions between 1.5 and 35 AU from the primary (A. Kraus 
and M. Ireland, private communication). The presence 
of ho t CO in the central 1 AU of the system (Salyk et al. 
|2007() and the high accretion rate, undiminished relative 
to the Taurus median, also argue against the presence 
of a massive close companion. A stellar companion is 
therefore an unlikely origin for the central clearing in the 
GM Aur system. 

Planet-Disk Interaction - Perhaps the most compelling 
mechanism for producing a transition disk is the dynamic 
clearing of material by a giant planet a few times the 
mass of Jupiter. The opening of gaps and holes in cir- 
cumstellar disks has long been pr edicted as a consequence 
of giant planet for mation (e.g. Lin fc Papaloizoufll986l : 
iBrvden et al.|[l999t ). Some simulations have shown that 
inner holes may in fact be a more common outcome than 
gaps as angular momentum transfer mediated by spiral 
density waves c an clear the inner disk faster than the vis- 
cous timescale (jVarniere et al.ll2006l : iLubow fc D'Angelo' 
I2OO6I ). The planet-induced clearing scenario w as consid- 
ered in detail for GM Aur bv lRice et all (|2003D and found 
to be globally consistent with the observed properties 
of the system (although their estimate of the inner hole 
radius is based on pie-Spitzer SED information). This 
mechanism naturally explains the diminished but per- 
sistent accretion rates and presence of small dust grains 



through two predictions of models of planet-disk inter- 
action: (1) filtration of dust grains according to size is 
expected at the inner disk edge, leading to a dominant 
popu lation of small grains in the inner disk (H ice et al.l 
'2006") ; and (2) a sustained reduction in accretion rate to 
~ 10% of that through the outer disk is predicted as the 
giant planet begins to intercept most of the accreting ma- 
terial l|Lubow fc D'Angelol2006D . These effects combined 
may also explain the enhanced gas-to-dust ratio in the in- 
ner disk. A planet-induced gap c ould also serve as a cat- 
alyst f or inside-out MRI clearing lChiang fc Murrav-Clavl 
(l20f)l . 

Given the observed 20 AU inner disk radius and the 
scenario of clearing via dynamical interaction with a gi- 
ant planet, it is possible to make a simple estimate of the 
distance of the planet from the star. The width of a gap 
opened by a plane t is approximately 2y/3 Roche radii 
(| Artvmowicd 1 1 9871 ) . and simulations show that the min- 
imum mas s necessary to open a gap is of order 1 Jupiter 
mass (e.g. lLin fc Pa oaloizou 1993: E dgar et al .''2007). If 
the outer edge of the planet-induced gap coincides with 
the 20 AU inner disk radius (with the portion of the disk 
interior to the planet cleared via spiral density waves or 
the MRI), then a companion between 1 and 40 times the 
mass of Jupiter would be located between 11 and 16 AU 
from the star. The influence of a planet carving out an 
inner cavity in the dust distribution is therefore a plau- 
sible scenario, bolstered by recent results demonstrating 
that a planet is responsible for dynamical sculpting of 
dust in the much older Fomalhaut system (jKalas et al.l 
[2008h . 

5.2. Evidence for a Warp? 

While the model comparison in fj4l above shows that 
CO emission from the disk is globally consistent with 
Keplerian rotation, the 11° difference in position angle 
between the continuum data and the two CO data sets 
is significant at the 5a level, and may indicate some 
kinematic deviation from pure Keplerian rotation in a 
single plane. Changes in position angle with physical 
scale are commonly interpreted as warps in the contex t 
of studies of galaxy dynamics (e.g. iRogstad et anil974f) : 
it may be that the change in position angle in the GM 
Aur disk indicates a kinematic warp. 

The possibility of a warp or oth er deviation fr om Ke- 
plerian rotation was discussed by Dutr ev et al.l |l998-), 
although their discussion was based on possible isophote 
twisting observed in integrated CO J=2— 1 contours. We 
observe no such isophote twisting in the integrated CO 
J =2-1 or J=3-2 emission presented here (Fig. |4]), al- 
though this determination may be influenced by the dif- 
fering baseline lengths and beam shapes in the respective 
interferometric data sets. Instead, we observe deviations 
from the expected position angle only in the rotation 
pattern of the resolved CO emission, which is reflected 
in the isovelocity contours of Fig. [H This position angle 
change does not appear to be related to the cloud con- 
tamination, as it is more clear in the less-contaminated 
CO J=3— 2 data set. In order to test whether the posi- 
tion angle of the true brightness distribution might have 
been altered by incomplete sampling of the data in the 
Fourier domain, we generated a model of the disk at a 
position angle of 64°, consistent with that measured in- 
dependently for the two continuum data sets. We then fit 
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the position angle by minimization as in above. 
With this method, after sampling with the response at 
the spatial frequencies in the CO J=3-2 data set, we 
recover the position angle to within less than a degree 
of the input model. This is to be expected, since the 
fitting procedure takes into account the interferom- 
eter response when fitting for the position angle. The 
position angle change is therefore robust independent of 
beam convolution effects. 

In order to cause a change in position angle on physical 
scales between those probed by the continuum (~ 30 AU) 
and the CO 200 AU), a warp would have to occur at 
a size scale of order 100 AU. The most natural explana- 
tions for the presence of a warp in a gas-rich circumstellar 
disk include flybys and perturbations by a planet or sub- 
stellar companion. A simple estimate of the timescale 
of flyby interactions is r = 1/ (Nirb'^a), where N is the 
number density of stars, b is the approach distance, and 
a is the velocity dispersion. Assuming typical values 
for Taurus, in cluding a stellar density of ~10pc~'^ (e.g. 
iGomez et al.|[l993i) and velo city dispersion of 0.2kms~^ 
( Kraus fc Hillenbrandl [20081 ) . the timescale for interac- 
tions at distances of ~1000 AU, sufficient to cause sig- 
nifica nt perturbations at Oort Cloud radii (jSchoU et al.l 
[l98l) . is of order 1 Gyr. Since the results of a one-time 
perturbation would likely damp in a few orbital periods 
(10^ yr at a distance of 100 AU), such an interaction is 
statistically unlikely. However, it should be noted that a 
recent interaction might have been capable of producing 
an extended feature like th e "blu e ribbon" observed in 
scattered light bv [Schneide r et al.l (|200 3). 

The influence of a massive planet or substellar com- 
panion has been investigated as the origin of warps 
observed in ga s -depl eted debris disks, inclu ding /? Pi c 
(jMouillet et al.l [l997l ) and HD 100546 (Quilie3 l200l) . 
However, there is a dearth of theoretical investigation 
into the plausibility of warps caused by planetary sys- 
tems in gas-rich disks more closely analogous to the GM 
Aur system. Since the warp in the GM Aur disk must 
occur between the Hill sphere of the putative planet and 
the ^^200 AU resolution of the CO line observations, it is 
plausible that the warp could be due to the gravitational 
influence of the same body responsible for evacuating 
the inner disk. A theoretical inquiry into this possibility 
would be useful, but is beyond the scope of this paper. 

6. CONCLUSIONS 

Spatially resolved observations in millimeter contin- 
uum emission, obtained using the SMA at 860 /zm and 
PdBI at 1.3 mm, reveal a sharp decrease in optical depth 
near the center of the GM Aur disk. Simple estimates of 
the extent of this region, based on the separation of peaks 
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in the continuum images and the position of the null in 
the visibility functions in Fig. (TJ are consistent w i th the 
inner hole radius of 24 AU derived bv lCalvet erall ()2005D 
using disk structure models to fit the SED. No significant 
azimuthal asymmetry is detected in the continuum emis- 
sion. 

Refined v e rsions of the SED-based model of 
ICalvet et al.l (|2005[ ) show that the data are very 
well reproduced by a disk model with an inner hole of 
radius 20 AU. This model overpredicts the broadband 
SED flux in the 20-40 ^m wavelength regime, but this 
region of the spectrum likely depends on the properties 
of the wall at the inner disk edge, which are poorly 
constrained by available data. 

CO emission in the J=3— 2 and J=2— 1 transitions 
conflrms the presence of a disk with kinematics consis- 
tent with Keplerian rotation about the central star, but 
at a position angle offset from the continuum by 
The morphology of the CO emission is broadly consistent 
with the SED model, but with a larger CO J=3--2/ J=2- 
1 line ratio than predicted for the SED model. This is 
a likely indication of additional gas heating relative to 
dust in the upper disk atmosphere. 

Given the observed properties of the GM Aur system, 
photoevaporation, grain growth, and binarity are un- 
likely physical mechanisms for inducing a sharp decrease 
in opacity or surface density at the disk center. The in- 
ner hole plausibly results from the dynamical influence 
of a planet on the disk material, with the inner disk pos- 
sibly cleared by spiral density waves or the MRI. While 
a recent flyby is statistically unlikely, warping induced 
by a planet could also explain the difference in position 
angle between the continuum and CO data sets. 
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